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Abstract 

We give the cumulative distribution functions, the expected values, and the mo- 
ments of weighted lattice polynomials when regarded as real functions of indepen- 
dent random variables. Since weighted lattice polynomial functions include ordinary 
lattice polynomial functions and, particularly, order statistics, our results encom- 
pass the corresponding formulas for these particular functions. We also provide an 
application to the reliability analysis of coherent systems. 
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1 Introduction 

The cumulative distribution functions (c.d.f.'s) and the moments of order 
statistics have been discovered and studied for many years (see for instance 
David and Nagaraja [4J). Their generalizations to lattice polynomial functions, 
which are nonsymmetric extensions of order statistics, were investigated very 
recently by Marichal [TU] for independent variables and then by Dukhovny [S] 
for dependent variables. 

Roughly speaking, an n-ary lattice polynomial is a well-formed expression 
involving n real variables Xi, . . . , x„ linked by the lattice operations A = min 
and V = max in an arbitrary combination of parentheses. In turn, such an 
expression naturally defines an n-ary lattice polynomial function. For instance, 

p{xi,X2, X3) = {xi A X2) V 0:3 
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is a 3-ary (ternary) lattice polynomial function. 

Lattice polynomial functions can be straightforwardly generalized by fixing 
certain variables as parameters, like in the 2-ary (binary) polynomial function 

p{xi, X2) = (c A Xi) V X2, 

where c is a real constant. Such "parameterized" lattice polynomial functions, 
called weighted lattice polynomial functions [7], play an important role in the 
areas of nonlinear aggregation and integration. Indeed, they include the whole 
class of discrete Sugeno integrals |12iil3j . which are very useful aggregation 
functions in many areas. More details about the Sugeno integrals and their 
applications can be found in the remarkable edited book [6J. 

In this paper we give closed-form formulas for the c.d.f. of any weighted lat- 
tice polynomial function in terms of the c.d.f. 's of its input variables. More 
precisely, considering a weighted lattice polynomial function p : M" — >■ M and 
n independent random variables Xi, . . . , X„, Xi {i = 1, . . . ,n) having c.d.f. 
Fi{x), we give formulas for the c.d.f. of Yp := p{Xi, . . . ,X„). We also yield 
formulas for the expected value E[5f(l^)], where g is any measurable function. 
The special cases g{x) = x, x^, [x — E(l^)]'', and e^^ give, respectively, the 
expected value, the raw moments, the central moments, and the moment- 
generating function of Yp. 

As the lattice polynomial functions are particular weighted lattice polyno- 
mial functions, we retrieve, as a corollary, the formulas of the c.d.f.'s and the 
moments of the lattice polynomial functions. 

This paper is organized as follows. In Section 2 we recall the basic material 
related to lattice polynomial functions and their weighted versions. In Section 
3 we provide the announced results. In Section 4 we investigate the particular 
case where the input random variables are uniformly distributed over the unit 
interval. Finally, in Section 5 we provide an application of our results to the 
reliability analysis of coherent systems. 



2 Weighted lattice polynomials 

In this section we recall some basic definitions and properties related to 
weighted lattice polynomial functions. More details and proofs can be found 
in [3. 

As we are concerned with weighted lattice polynomial functions of random 
variables, we do not consider weighted lattice polynomial functions on a gen- 
eral lattice, but simply on a closed interval L := [a, b] of the extended real 
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number system M := M U {— c)o, +00}. Clearly, such an interval is a bounded 
distributive lattice, with a and h as bottom and top elements. The lattice 
operations A and V then represent the minimum and maximum operations, 
respectively. To simplify the notation, we also set [n] := {l,...,n} for any 
integer n ^ 1. 

Let us first recall the definition of lattice polynomial functions (with real 
variables); see for instance Birkhoff [2, §11.5]. 

Definition 1 The class of lattice polynomial functions from to L is de- 
fined as follows: 

(1) For any k G [n], the projection {xi, . . . , x„) ^ Xk is a lattice polynomial 
function from to L. 

(2) If p and q are lattice polynomial functions from L" to L, then p A q and 
pV q are lattice polynomial functions from to L. 

(3) Every lattice polynomial function from to L is constructed by finitely 
many applications of the rules (1) and (2). 

As mentioned in the introduction, weighted lattice polynomial functions gener- 
alize lattice polynomial functions by considering both variables and constants. 
We thus have the following definition. 

Definition 2 The class of weighted lattice polynomial functions from L" to 
L is defined as follows: 

(1) For any k G [n] and any c & L, the projection {xi, . . . ,Xn) ^ Xk and 
the constant function 1— c are weighted lattice polynomial 
functions from to L. 

(2) If p and q are weighted lattice polynomial functions from to L, then 
p A q and pM q are weighted lattice polynomial functions from to L. 

(3) Every weighted lattice polynomial function from to L is constructed 
hy finitely many applications of the rules (1) and (2). 

Because L is a distributive lattice, any weighted lattice polynomial function 
can be written in disjunctive and conjunctive forms as follows. 

Proposition 3 Let p : L"- —>■ L be any weighted lattice polynomial function. 
Then there are set functions a : 2^^^ L and (3 : 2^^^ — »• L such that 



p(x)= V 


a{S) A f\xi 


= A 


I3{S) V V ^^ 


SC[n] 




SC[n] 


'- i&S 



Proposition [3] naturally includes the classical lattice polynomial functions. To 
see this, it suffices to consider nonconstant set functions a : 2'"] {a, h] and 
(5 : 2W ^ {a, 6}, with a(0) = a and /?(0) = h. 
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The set functions a and (5 that disjunctively and conjunctively define the 
polynomial function p in Proposition [3] are not unique. For example, we have 



X\ V [x\ A X2) = X\ = X\ f\ {x\ V X2). 



However, it can be shown that, from among all the possible set functions 
that disjunctively define a given weighted lattice polynomial function, only 
one is nondecreasing. Similarly, from among all the possible set functions that 
conjunctively define a given weighted lattice polynomial function, only one is 
nonincreasing. These particular set functions are given by 

a{S) = 'p{&s) and /5(S') = p(e[„]\5), 

where, for any C [n], 65 denotes the characteristic vector of S in {a, 6}", i.e., 
the n-dimensional vector whose ith component is 6, if i G S", and a, otherwise. 
Thus, any n-ary weighted lattice polynomial function can always be written 
as 



5CM 



A P(eN\5)vV 



The best known instances of weighted lattice polynomial functions are given by 
the discrete Sugeno integrals, which consist of a nonlinear discrete integration 
with respect to a fuzzy measure. 

Definition 4 An L-valued fuzzy measure on [n] is a nondecreasing set func- 
tion /i : 2["] L such that /i(0) = a and fi{[n]) = b. 

The Sugeno integrals can be presented in various equivalent forms. The next 
definition introduce them in one of their simplest forms (see Sugeno ^2\)- 



Definition 5 Let fi be an L-valued fuzzy measure on [n] . The Sugeno integral 
of a function x : [n] — L with respect to fi is defined by 



SC[n] 



Thus, any function / : ^ L is an ra-ary Sugeno integral if and only if it is 
a weighted lattice polynomial function fulfilling /(e^) = a and /(e[„]) = b. 



3 Cumulative distribution functions and moments 



Consider n independent random variables Xi,...,X„, Xj {i G [n]) having 
c.d.f. Fi{x), and set Yp := p{Xi, . . . ,X„), where p : L" — L is any weighted 



4 



lattice polynomial function. Let H : W ^ {0, 1} be the Heaviside step function, 
defined by H{x) = 1, if a; ^ 0, and 0, otherwise. For any c G M, we also define 
the function Hc{x) := H{x — c). 

The c.d.f. of Yp is given in the next theorem. 

Theorem 6 Let p : ^ L he a weighted lattice polynomial Junction. Then, 
the c.d.f. ofYp is given by each of the following formulas: 



= [l-//,(e,)(y)] n F,{y) nil - F,{y)], 

SC[n] i&[n]\S if^S 

SC[n] ieS ie[n]\S 



(1) 

(2) 



Proof. Starting from the disjunctive form of p, we have 



Fp{y) = l-FT 
= 1-Pr 



V [pies)A/\X, 

SC[n] 



> y 



3S C [n] such that y < p{es) and y < /\ 

ies 



Consider the following events: 



A:-- 
B:-- 



3S C [n] such that y < p{es) and y < f\Xi 



ies 



3S C [n] such that y < p{es) and \/ X^ ^ y < /\ X^ 

iG[n]\S i€S 



Event B implies event A trivially. Conversely, noting that p is nondecreasing 
in each variable and replacing S with a superset S' I) S, if necessary, we 
readily see that event A implies event B. 



Since the events 
we have 



Vie[n]\5^i ^y < Kes^i {S Q [n]) are mutually exclusive. 



Fp{y) = l- E [y < Pi^s) 

SC[n] 



Ft 



V X,^y< /\Xi 
ie[n]\s ieS 



which, using independence, proves the first formula. The second one can be 
proved similarly by starting from the conjunctive form of p. □ 
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The expressions of Fp{y), given in Theorem [6l are closely related to the follow- 
ing concept of multilinear extension of a set function, which was introduced 
by Owen ^llj in game theory. 

Definition 7 The multilinear extension of a set function v : 2t"'l ^ M is the 
function $^ : [0, l]** —>■ M defined by 

$,(x):= Yl viS)l[x, n (l-a^O- 

SC[n] ieS ie[n]\S 



Using this concept, we can immediately rewrite ([T]) and ([2]) as 

F,{y) = 1 - <|.„^J1 - F,{y), K{y)], 
F,(y) = $,,jFi(|/),...,F„(|/)], 

where, for any fixed y E M, the (nondecreasing) set functions Vp^y : 2'" 
{0, 1} and v^y : 2^ ^ {0, 1} are defined by 

Vp,y{S) := 1 - i^p(es)(l/) and v^y^S) := Hp(^^^^^^^){y). 



Owen [TT] showed that the function being a multilinear polynomial, has 
the form 

$,(x) = J2 ^-(s) E^i, (3) 

SC[n] i€S 

where the set function ruy : 2'"] — M, called the Mobius transform of f , is 
defined as 

^.(^)= E(-i)'"""^'^(^)- (4) 

TCS 

Using this polynomial form of we can immediately derive two further 
formulas for Fp{y), namely 



Fp{y) = 1 - E ^^.yiS) nil - FM], (5) 

SCln] i&S 

Fp{y) = Y: my*JS) l[F^iy). (6) 

5C[n] ieS 



Formulas ([I])-([2]) and ([S])-© thus provide four equivalent expressions for 
Fp{y). As particular cases, we retrieve the c.d.f. of any lattice polynomial 
function. For example, using formula ([1]) leads to the following corollary (see 

m)- 
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Corollary 8 Let p : ^ L be a lattice polynomial function. Then, the c.d.f. 
of Yp is given by 

FM = ^- E n F^iy) Ul^ - F^iy)]- 

SC[n\ i€ln]\S ieS 
p{es)=b 

Let us now consider the expected value E[(yf(yp)], where : R — > R is any mea- 
surable function. For instance, when g{x) = x^, we obtain the raw moments 
ofYp. 

By definition, we simply have 

/OO POO 
g{y)dFp{y) = - g (y) d[l - Fp{y)]. 
-OO J — OO 

Using integration by parts, we can derive alternative expressions of £[(7(1^)]. 
We then have the following immediate result. 

Proposition 9 Let p : L be any weighted lattice polynomial function 

and let g : W ^ K. be any measurable function of bounded variation. 

(1) If\\vsiy^^g{y)[l - Fiiy)] = for all i G [n], then 

/OO 
[1-Fp{y)]dg{y). 
-00 

(2) If\m\y^_^g{y)F.i{y) = for alii E [n], then 

/OO 
Fp{y)My)- 
-00 

Clearly, combining Proposition [H] with formulas and ([S])-® immedi- 

ately leads to various explicit expressions of £[(7(1^)]. For instance, if 

hmgiy)[l-F,iy)] = 0, G [n], 

then 




It is noteworthy that Eq. ([7]) can also be established without the knowledge 
of the c.d.f. of Yp. As the proof is very informative, we give it in the appendix. 

Remark 10 We can also retrieve the c.d.f. of Yp directly from formula 
Indeed, rewriting (0j with the function g{y) = H{z — y), we immediately obtain 
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= l+r E [l-^P(es)(l/)] n F,{y)J{[l-F,{y)]dH{z-y) 



5C[n 

and hence we retrieve H 



i&[n]\S 



Remark 11 Suppose that some variables Xi are constants, say, Xk = Ck for 
all k & K for a given K ^ [n]. Then, the weighted lattice polynomial function 
p reduces to an {n — \K\)-ary weighted lattice polynomial function and it is 
easy to see that 



V p{es)^^X, 

SC[nl ieS 



V a^{S) A/\X, 

SC[n]\K i£S 



where 



a, 



K/ 



TCK 



{S) := V Pi^sur) A A c, 



{SC[n]\K). 



Thus, the conditional expectation Fj[g{Yp) \ X^ = Ck^k ^ K] can be immedi- 
ately calculated by Proposition^^ 



4 The case of uniformly distributed variables on the unit interval 



We now examine the case where the random variables Xi, . . . , X„ are uni- 
formly distributed on [0, 1]. We also assume L = [0, 1]. 

Recall that the incomplete beta function is defined, for any u,v > 0, by 

BJu,v) := r f'-^l-ty-^dt (zeR), 
Jo 

and the beta function is defined, for any u,v > 0, by B{u, v) := Bi{u, v). 

According to Eq. for any weighted lattice polynomial function p : [0, 1]" — 
[0, 1] and any measurable function : [0, 1] — R of bounded variation, we have 

ngiY,)] = g{o) + E / - yf^ My)- 

If, furthermore, g{y) = ^y^, with r G N \ {0}, then 

;E[r;]= E B,^,,)in-\S\+r,\S\ + l). (8) 

SC[n] 
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Remark 12 Considering Eq. with p{es) = z, if S = [n], and 0, other- 
wise, we obtain 



BJr,n+l] 



and hence the following identity 



1 

r J [0,1]" 



r J[o,i]'^ 



SC[n] 



V pMaA 



Xi 



i€S 



z A /\ Xi 



1 



dx 



P(e5)A/\ Xi 

i€S 



n-\S\+r 



which shows that computing the raw moments of any weighted lattice polyno- 
mial reduces to computing the raw moments of bounded minima. 

Let us now examine the case of the Sugeno integrals. As these integrals are usu- 
ally considered over the domain [0, 1]", we naturally calculate their expected 
values when their input variables are uniformly distributed over [0, 1]". Since 
any Sugeno integral is a particular weighted lattice polynomial, by Eq. ([8]), its 
expected value then writes 



/ 5^(x)dx= E B^(^s){n-\S\ + l,\S\ + l) 

•^10.11" SQn] 

V / x"-l^l(l-a;)l^ldx 



SC[n] 

SC[n] fi=0 V ^ / 



n-\S\+i + l 



Surprisingly, this expression is very close to that of the expected value of the 
Choquet integral with respect to the same fuzzy measure. 

Let us recall the definition of the Choquet integrals [3j. Just as for the Sugeno 
integrals, the Choquet integrals can be expressed in various equivalent forms. 
We present them in one of their simplest forms (see for instance [8]). 

Definition 13 Let ^ be a [0,l]-valued fuzzy measure on [n]. The Choquet 
integral of a function x : [n] ^ [0, 1] with respect to /i is defined by 



E ^t^(^) A 

SOn] ie5 



where m,, : 2'"' 



zs 



the Mobius transform (cf. ^) of ^. 



For comparison purposes, we recall the expected value of (see for instance 
0): 
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X 



[0,1] 



C^(x) dx 



^ fi{S)B{n-\S\ + l,\S\ + l) 




5 Application to reliability theory 

In this final section we show how the results derived here can be applied to the 
reliability analysis of coherent systems. For a reference on reliability theory, 
see for instance Barlow and Proschan [T]. 

Consider a system made up of n independent components, each component Ci 
{i G [n]) having a lifetime Xi and a reliability rj(t) := Pr[Xj > t] at time t > 0. 
Additional components, with constant lifetimes, may also be considered. 

We assume that, when components are connected in series, the lifetime of 
the subsystem they form is simply given by the minimum of the component 
lifetimes. Similarly, for a parallel connection, the subsystem lifetime is the 
maximum of the component lifetimes. 

It follows immediately that, for a system mixing series and parallel connec- 
tions, the system lifetime is given by a weighted lattice polynomial function 



of the component lifetimes. Our results then provide explicit formulas for the 
c.d.f., the expected value, and the moments of the system lifetime. 

For example, the system reliability at time t > is given by 



Moreover, for any measurable function g : [0,oo] ^ M of bounded variation 
and such that limy^oo 9 {y)ri{y) = for all i G [n], we have, by Proposition [9l 



As an example, the following proposition yields the mean time to failure E[Yp] 
in the special case of the exponential reliability model. 



i; = p(Ai,...,x„) 



Rp{t) := Pi[Yp >t] = $.,,Jri(t), . . . ,r„(t)]. 



(9) 
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Proposition 14 If ri{t) = e for alii G [n], we have 

SC[n] TCS '^WJ 

where A (5) := Eies -hi- 



proof. Using Q) and then ([SD^dll), we obtain 

Rp{t)= E^US)e~'^'^^= E E(-l)'"'~'^'-pXr)e-^(^)* 

5Cfnl SCfnlTCS 



and lience 

5c[n] res 



E 



POO pp(eT) 

[Y,]= / R,{t)dt= E E(-l)'"""^' / e-^(^)Mt. □ 



6 Conclusion 



We have extended the c.d.f.'s and moments of lattice polynomial functions 
to the weighted case. At first glance, this extension may appear as a simple 
exercise. However, it led us to nontrivial formulas, which can be directly ap- 
plied to qualitative aggregation functions such as the Sugeno integrals and 
their particular cases: the weighted minima, the weighted maxima, and their 
ordered versions. 



Appendix: Alternative proof of Eq. ([7]) 



In this appendix we present a proof of Eq. ([7]) without using the c.d.f. of 
Yp. The main idea of this proof is based on the fact that the expected value 
E[(7(y^)] can be expressed as an n-dimensional integral, namely 

/oo roo 
•■■/ g[p{^)]dF,{xi)---dFn{xn). (10) 
-oo J — oo 



At first glance, the evaluation of this expression requires a difficult or in- 
tractable integration. However, as we will now see, any weighted lattice poly- 
nomial function fulfills a remarkable decomposition formula, which will enable 
us to calculate E[5f(l^)] in a straightforward manner. 
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Given a weighted lattice polynomial function p : ^ L and an index A; G [n], 
define the weighted lattice polynomial functions pf.-.L'^—^L and p\: ^ L 

as 



Pfc(x) ■.= p{Xi, . . . ,Xk-l,a,Xk+l, . . .,Xn), 
pI{x) ■.= p{Xi, . . .,Xk-l,b, Xk+l, . . .,Xn). 

Then, it can be shown [7j that 



p(x) = median 



Pfc(x),Xfc,p^(x) 



UeL^). fill 



This decomposition formula expresses that, for any index k, the variable Xk can 
be totally isolated in p(x) by means of a median calculated over the variable 
Xk and the two functions p^ and p^., which are independent of Xk- 

This interesting property leads to the following lemma. 

Lemma 15 Let p : ^ L he any weighted lattice polynomial function, let 
G [n], and let g : W ^ be any measurable function of bounded variation 
and such that g{—oo) is finite and liiriy g{y)[i — Fk{y)] = 0. Then 



E[giY^] = E[ga{Y,.)] + E[g,iY^.)], (12) 



where 



ga{x):= r Fk{t)dgit), 

J —oo 

g,{x):=g{x)- f Fk{t)dg(t). 

J —CO 



Proof. Let k G [n] and fix Xj for all j ^ k. Assume u := ^^(x) and v := Pfc(x) 
are finite. The other cases can be dealt with similarly. Then, we have u ^ v 
and, by (fTTl) . 

p(x) = median[M, Xk, v]. 
Hence, for any measurable function (7 : M ^ M of bounded variation, we have 

/oo ru rv roc 

g[p{x)]dFk{xk)= g{u)dFk{xk)+ g{xk) dFk{xk) + g{v)dFk{xk) 
-00 J —00 Ju Jv 

= g{v)- f Fk{t)dg{t) 
= 9aiu) + 9biv)- 

Finally, integrating with respect to the other variables Xj, j 7^ k, and then 
using (ITU]) , we get the result. □ 
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Observe that formula (fT2!) . when considered for every index k and every func- 
tion g, completely determines the expected value of g(Yp). Indeed, repeated 
applications of that formula will eventually lead to integration of transformed 
weighted lattice polynomial functions all of whose variables are set to either 
a or b. 



Proof of Eq. ([7]). For any fixed S C [n], define recursively the sequence 
{gsYj^^Q of functions g^ : R ^ M as g^ := g and, for k ^ 1, 




Fk{t)dgl~\t), ifk^S, 
- r F,{t)dgl-\t), ifkeS. 

Repeated applications of Lemma [15] eventually lead to 

SC[n] 



Let us now show that 



gUz) = 9si-oo) + r n Fk{x) Ul^ - F,{x)] dg{x) 



For any 5* C [n] and any k E [n], we have 

Fu{x)dg^s~\x), iikiS 



^9s{x) 

and hence 



[l-Fk{x)]dgt\x), likes, 
d^s(^)= n F,{x)Xl[l-F,{x)]dg{x), 

k&[n]\S k£S 



which proves the result since gg{—oo) = g{—oo), if S = N, and 0, other- 
wise. □ 
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